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been  determined.  A  new  method  for  calculating  electronic 
autoionization  widths  using  high  Rydberg  state  wave  functions  to 
represent  the  inner  part  of  the  free  electron  wave  function  has 
b^n- used  and^  tested  -on  the_NO  molecule  wher^we  have  found  that 
the  widths  can  be  calculated  to  arT'^pected^^accuracy  of  about 
15%.  The  widths  have  been  used  to  determine  dissociative 
recombination  cross  sections  as  a  function  of  electron  energy. 
Large  windows  have  been  discovered  in  the  cross  sections  from 
excited  vibrational  levels.  The  windows,  at  which  the  cross 
sections  drop  precipitously,  are  due  to  the  overlap  of  the  peak 
in  the  continuum  vibrational  wave  f unctions the  number  of 
windows  in  the  cross  section  from  a  particular  ion  vibrational 
level  is  given  by  the  number  of  ion  vibrational  wave  function 
nodes  lying  to  the  small  internuclear  distance  side  of  the 
intersection  point  of  the  repulsive  curve  with  the  vibrational 
level  undergoing  recombination.  A  single  dissociative  route  has 
been  identified  as  the  main  channel  leading  to  excited  atoms 
while  three  routes  have  been  shown  to  be  important  for  the 
production  of  atoms.  The  production  of  S  atoms  is  very 
sensitive  to  vibrational  excitation.  The  rate  fr<xn  the  third 
vibrational  level  is  about  80  times  larger  than  the  rate  from  the 
lowest  vibrational  level.  For  atoms  the  irate  from  v*0  is  the 
largest  and  is  about  a  factor  of  two  larger  than  the  rate  from 


mnagsmii 


Jliuwwuviw 


nuvwi 


kS*^' 


»v»' 


MOLECULAR  SOURCES  OF 
IONOSPHERIC  HOLES 


AFOSR-84-0109 


.TH-  8  8^0064 


Final  Report 

For  the  period  1  April  1984  through  30  September  1987 


Principal  Investigator 
Dr.  Steven  L.  Guberman 


November  1987 


Prepared  for 

Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  D.  C.  20332 


Smithsonian  Institution 
Astrophysical  Observatory 
Cambridge,  Massachusetts  02138 


The  Smithsonian  Astrophysical  Observatory 
is  a  member  of  the 

Harvard-Smithsonian  Center  for  Astrophysics! 


The  AFOSR  Technical  Officer  for  this  grant  is  Major  Larry  Davis, 
Directorate  of  Chemical  and  Atmospheric  Sciences,  AFOSR/NC, 

Bolling  AFB,  D.  C.  20332 


88  2 


A 

.-V 


18  8 


k 


TABLE  OF  CONTENTS 


I.  Introduction 

1 

II.  Wave  Functions  and 
Potential  Curves 

2 

III.  Electronic  Widths 

5 

IV.  Cross  Sections 

10 

V.  References 

13 

VI.  Publications 

15 

VII.  Presentations 

16 

VIII.  Figure 

18 

IX.  Appendix 

21 

Accession  For 

NTIS 

GR.'-il 

m 

DTI-' 

TAH 

n 

'  .It*  ,;i 

.'.'d 

□ 

i  4  '  '  ! 

*  '  ‘a  _  1. 

a - - - 

A  , 

L- 

‘  ^  Ul  1  OL'/ 

AvMi'lblllty  ro;1C;3 

A  ve.l  1 

11'  j/ or 

.IS. 

S'p^c 

'.'ll 

r' 

FINAL  TECHNICAL  REPORT 


AFOSR-8i4-O109 

MOLECULAR  SOURCES  OF  IONOSPHERIC  HOLES 
For  the  period  April  1,  1984  -  September  31 f  1987 


I.  Introduction 


The  dissociative  recombination  (DR)  of  molecular  ions  with  electrons  is  an 
important  process  determining  the  radiative  characteristics  and  overall 
energy  balance  of  both  the  quiescent  and  perturbed  ionosphere.  The  large 
scale  ionospheric  electron  depletion  or  hole  that  results  after  the  launch 
of  large  rockets  has  been  attributed  to  DR.*  The  modelling  of  the  detailed 
role  played  by  DR  in  the  Earth's  atmosphere  has  been  hampered  by  the  lack  of 
reliable  rate  constants.  During  the  past  3  1/2  years  we  have  performed  large 
scale  ab  initio  calculations  of  the  rate  constants  for  the  DR  of  Oj  with  an 
electron: 


Ot  +  e~  -«■  0  +  0.  (1) 

The  research  has  combined  state  of  the  art  techniques  in  both  bound  state 
calculations  and  scattering  theory.  The  results  provide  a  detailed 
description  of  the  variation  of  the  DR  rate  constants  for  producing  excited 
states  of  0  with  ion  vibrational  excitation  and  electron  temperature.  In 
addition  to  being  an  application  to  Oj  DR  the  research  has  involved  the 


development  of  new  techniques  for  the  accurate  calculation  of  molecular 
electronic  autoionization  widths  and  potential  energy  curves.  Results 


obtained  from  these  new  techniques  are  summarized  in  Sections  II  and  III 
below.  Section  IV  presents  the  cross  section  calculations  and  a  discussion 
of  the  large  windows  which  have  been  discovered  in  DR  cross  sections  from 
excited  ion  vibrational  levels.  Results  for  the  generation  of  both  0(^S)  and 
0{*D)  atoms  from  DR  are  presented  and  compared  to  satellite  observations  in 
Section  IV  and  the  Appendix. 

II.  Wave  Functions  and  Potential  Curves 


Potential  curves  determined  previously*'*  have  been  improved  by  using 
considerably  larger  basis  sets  and  configuration  interaction  (Cl)  wave 
functions  having  a  factor  of  100  more  terms  than  in  the  earlier  studies. 
Briefly,  the  Cl  wave  function  is  composed  of  configurations,  each  of  which 
describes  a  possible  distribution  of  the  electrons  in  the  molecular  orbitals 
for  that  electronic  state.  For  example,  for  the  ground  state  of  Oj,  the  most 
important  configuration  near  the  equilibrium  internuclear  separation  is: 

(10g)^(1  Ou)*(20g)*(20y)*(30g)*(lTtu)‘'  (iTTg)*  (2) 

The  la  and  2a  orbitals  are  mostly  0  atom  Is  and  23  orbitals,  respectively. 
The  3og  orbital  is  a  bonding  combination  consisting  mostly  of  the  0  atom  2p 
orbitals  lying  along  the  bonding  axis.  The  and  Wg  orbitals  are  out  of 
plane  bonding  and  antibonding  orbitals  respectively,  consisting  mostly  of 
the  0  2p  orbitals  that  are  perpendicular  to  the  bonding  axis.  In  the 
techniques  used  here,  the  orbitals  in  (2)  are  optimized  for  the  electronic 
state  of  interest  (in  this  example  the  ground  state  of  0^)  in  a  large 
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calculation  which  includes  all  the  configurations  of  the  proper  total 
symmetry)  which  can  be  formed  by  distributing  the  electrons  in  all  the 
valence  orbitals.  The  latter  set  would  include  the  3ou  oi'bital  which  is  not 
Included  in  (2).  In  general,  the  1o  orbitals,  which  do  not  participate  in 
the  bonding  are  constrained  to  remain  doubly  occupied.  For  the  ground  state 
of  Oj  this  procedure  generates  14  terms  or  configurations.  Each  of  the 
orbitals,  <J>j,  including  those  shown  in  (2),  are  expanded  in  a  large  set  of 
Gaussian  basis  functions,  xi« 


♦j  ■  ^  ‘^J.i  Xi  ^3) 

1 

Both  the  coefficient  in  front  of  each  configuration  and  the  cj^i  are 
optimized  in  multiconfiguration  self  consistent  field  (MCSCF)  or  in  Complete 
Active  Space  Self  Consistent  Field  (CASSCF)  calculations.**  In  general,  many 
of  the  Xi  In  (3)  consist  of  contractions  of  individual  or  primitive 
Gausslans®  whose  weighting  coefficients  have  been  determined  in  atomic 
calculations.®  Since  the  orbitals  are  expanded  in  large  basis  sets,  Xii  the 
number  of  orbitals  used  in  the  MCSCF  or  CASSCF  is  only  a  subset  of  the  full 
number  of  orbitals  that  can  be  formed  from  the  basis  set.  The  remaining 
orbitals  are  called  virtual  orbitals  and  are  employed  in  the  Cl  calculations 
to  obtain  even  more  accurate  results  than  those  from  the  MCSCF  or  CASSCF 
procedures.  The  Cl  wave  functions  are  constructed  by  exciting  electrons  out 
of  a  small  set  of  reference  configurations  into  the  virtual  orbitals.  The 
reference  configurations  usually  consist  of  all  the  conf Igurations  needed  to 
properly  dissociate  that  particular  state  to  the  proper  separated  atoms 
limit.  Taking  at  most  one  or  two  electrons  from  a  reference  configuration 
and  promoting  them  to  the  virtual  orbitals  leads  to  a  singles  and  doubles  Cl 
denoted  CISD.  The  calculation  of  these  large  wave  functions  often  having 
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more  than  100,000  terms  has  been  made  possible  in  recent  years  by  the 
development  of  direct  Cl  techniques  for  solving  for  the  mixing  coefficients 
in  front  of  each  configuration.’  The  calculations  used  here  employ  these 
large  scale  wave  functions  and  direct  Cl  techniques.  Calculations  on  several 
states  of  Oj  and  N*  provide  examples  of  the  sort  of  accuracy  which  can  be 
expected  with  these  techniques.  Calculations  on  the  ground  state  of  O2  using 
a  basis  set  consisting  of  6s,3p,2d  and  If  type  contracted  Gaussians  leads  to 
a  228,036  term  wave  function.  The  calculated  (experimental*)  results  are 
1565. '17(1580. 19)  for  ue.  10 .93(  1 1  .98cm~M  for  ueXe.  and  2. 2999(2. 28l9ao)  for 
Rg.  For  Oj,  a  series  of  vertical  excitation  energies  have  been  calculated. 
The  calculated  (experimentally  derived)  results  are  5. 8982(5. 90eV)  for  *1^, 
6. 1318(6. lleV)  for  6. 2766(6. 19eV)  for  ’Z*,  and  8.66(8.6l±0.05eV*)  for 
’ZyT  All  the  experimental  energies  are  derived*®  from  Ref. 8  except  for  “ZyT 
These  calculated  vertical  excitation  energies  deviate  by  no  more  than  0.09eV 
from  experiment  and  indicate  that  the  placement  of  the  calculated  repulsive 
curves  relative  to  the  ion  should  be  highly  accurate.  For  the  ground  state 
of  Nj ,  14  reference  configurations  (needed  to  properly  dissociate  the  ground 
state  to  "S  atoms)  are  used  in  the  Cl  with  a  basis  set  of  the  same  size 
leading  to  a  Cl  wave  function  having  176,536  terms.  The  calculated 
(experimental**)  spectroscopic  constants  are  2340. 57(2356. 56cm“M  for  wg, 

14. 27(14. 3244cm“1 )  for  ugXg,  and  2.0856(2 .O74347ao)  for  the  equilibrium 
internuclear  distance,  Rg. 

The  techniques  described  here  have  been  applied  to  the  Important  DR 
routes  for  the  low  ion  vibrational  levels.  The  results'*  indicate  that  a 
state  of  * Zy  symmetry  which  crosses  the  ion  between  the  large  R  turning 
points  of  v-1  and  v-2  is  the  only  state  which  leads  to  0(*S)  atoms  from  the 
lowest  10  vibrational  levels.  On  the  other  hand,  potential  curve 
considerations  alone  indicate  that  there  are  five  states  that  can  contribute 


to  0(‘D)  formation  from  the  lowest  ion  vibrational  levels:  1 ‘llg,  I’lJ, 

and  The  results  are  described  in  more  detail  in  the  Appendix 

and  Section  IV. 

III.  Electronic  Widths 


The  calculation  of  DR  cross  sections  requires  the  calculation  of  the 
auto ionization  lifetimes  or  widths  of  the  molecular  states  which  cross  into 
the  electron-ion  continuum.  Techniques  for  calculating  the  widths  have  been 
developed  with  AFOSR  support.  The  calculation  of  molecular  electronic 
capture  widths  has  historically  been  a  difficult  problem  for  quantum 
chemists.  The  width,  f,  is  a  matrix  element  of  the  Hamiltonian  operator,  H, 
given  by  Fermi's  golden  rule: 


r(R)  =  2iTp|<yf>|H|'r5>|2 


(4) 


where  is  the  wave  function  of  the  properly  antisymmetrized  product  of  a 
"free"  electron  in  a  coulomb  like  orbital  and  the  wave  function  for  the  ion. 
Vf  is  the  electronic  wave  function  for  the  neutral  state  which  describes  the 
products  of  reaction  (1)  or  the  entrance  channel  of  reaction  (4),  R  is  the 
internuclear  distance,  and  p  is  the  density  of  states.  In  order  to  do  the 
integral  in  (4)  both  the  coulomb  orbitals  and  the  bound  orbitals  must  be 
expanded  about  a  common  origin.  With  the  coulomb  orbital  placed  at  the 
molecular  midpoint,  the  bound  orbitals,  which  have  most  of  their  amplitude 
near  the  nuclei,  must  instead  be  represented  at  the  molecular  midpoint  by  a 
very  large  slowly  convergent  basis  function  expansion.  Similar  difficulties 
arise  in  solving  for  the  optimum  coulomb  orbital  in  the  field  of  the  bound 
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electrons.  Under  AFOSR  support,  these  problems  have  been  overcome  by  using 
Rydberg  orbitals  with  high  principal  quantum  numbers,  n,  to  describe  the 
coulomb  orbital.  Except  for  the  coulomb  orbital,  all  the  orbitals  in  (^)  are 
bound  and  have  large  amplitude  only  near  the  molecule.  As  a  result,  it  is 
only  important  to  know  the  amplitude  of  the  coulomb  orbital  near  the 
molecule.  The  inner  part  of  coulomb  orbitals  are  well  known  to  be  very 
similar  to  the  inner  part  of  high  Rydberg  orbitals.*'*  With  the  coulomb 
orbital  represented  as  a  high  Rydberg  orbital  all  the  orbitals  in  (4)  are 
then  bound  and  are  represented  in  terms  of  Gaussian  basis  functions.  Large 
slowly  converging  expansions  are  not  needed.  Since  Gauss ians  have  the 
property  that  the  product  of  two  Gaussians  on  different  centers  is  just  a 
third  Gaussian  located  along  the  line  connecting  the  first  two,  the  matrix 
element  in  (^)  can  be  easily  evaluated. 

The  Rydberg  orbitals  for  the  n  states  of  NO  have  been  represented  by  a 
set  of  18  diffuse  Gaussians  basis  functions  placed  at  the  molecular 
midpoint.  The  Rydberg  orbitals  are  determined  in  Improved  Virtual  Orbital** 
calculations  in  the  field  of  the  ground  state  of  the  ion.  Rydberg  states 
up  to  n-8  can  be  accurately  represented.  Using  the  orbitals  determined  in 
this  procedure,  configuration  interaction  (Cl)  wave  functions  for  the  ion 
and  neutral  states  are  then  constructed.  The  Rydberg  states  are  further 
optimized  in  the  field  of  an  optical  potential  due  to  the  dissociating 
states  and  the  matrix  element  in  (4)  is  calculated  directly  from  the 
optimized  wave  functions.  The  procedure  is  outlined  below. 

For  n*>>*.,  where  il  is  the  angular  momentum  quantum  number  it  can  be 
shown**  that  the  matrix  element  in  (4)  is  given  by 


|<'l'f  |H|4'f>  1^  =  ( 1 /n*3)kexp(-c/n*) 


(5) 


*«* 


»  .  *  Mil  .yi  A 


^vWTWjwcwrar^  r 


vrtiere  k  and  c  are  constants  determined  from  the  calculation  of  the  matrix 
elements  on  the  left  side  of  (5)  for  the  highest  optimized  Rydberg  states. 
The  high  calculated  Rydberg  states  are  nearly  hydrogenic  and  the  energies 
can  be  represented  as  E=*-(  1 /2(n* )^)  where  n*=n-6  and  6  is  the  quantum 
defect.  The  density  of  states  can  be  written  as 

p  =  1/(E(n*'-1/2)  -  E(n*-H/2))  (6) 

which  can  be  closely  represented  as 

p  =  n*3exp(-1/(2n*2)).  (7) 

Substituting  (7)  and  (5)  into  (*4)  gives 

r  »  27rkexp( -(c/n*)-( . 5/n*^) ) .  (8) 

Taking  the  limit  of  n*  -»  ®  gives  the  threshold  capture  width, 

r  =  2TTk.  (9) 


The  Rydberg  wave  functions  have  been  refined  further  by  optimizing  them  in 
the  field  of  an  optical  potential  due  to  the  'if  states.  Feshbach  projection 
operators’^  have  been  defined  which  partition  the  total  space  such  that  y  = 
Py  +  Qy  where  Py  contains  the  Rydberg  states  to  be  optimized  and  Qy  are  the 
'if  states  described  above.  In  the  usual  Feshbach  projection  operator 
formalism,  an  eigenvalue  equation  is  derived  for  Py  that  is  difficult  to 
solve  because  it  contains  an  energy-dependent  optical  potential.  However,  a 
new  technique’®,  based  on  earlier  advances,'®  has  been  developed  in  which  Py 
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can  easily  be  determined  by  solving  the  usual  Cl  problem.  Writing  the  total 
wave  function  as  'i'  =  PY  +  QH*,  we  can  write  the  Schroedinger  equation  in 
matrix  form  as 


(10) 


where  Hpp  =•  PHP,  Hgp  =  QHP,  etc. 

Multiplying  the  matrices  in  (10)  leads  to: 


MppP*!*  +  HpqQf  =  EPy 

and 

HqpPf  +  HQQQy  =  EQH'. 

From  Eq.(12)  we  have, 


(11  ) 


(12) 


QH-  =  HQpPf/(E  -  Hqq). 


(13) 


Substituting  Eq.(13)  into  Eq.(11)  leads  to  a  matrix  optical  potential  for 

py, 


(Hpp  +  HpqHqp/(E  -  Hqq))  py  =  EPy.  (14) 

Since  E  is  on  both  sides  of  Eq.(l4),  it  is  difficult  to  solve  directly  for 
Py.  However  Py  can  be  easily  determined  by  simply  diagonalizing  the  full  H 
matrix  in  Eq.(IO)  and  retaining  only  the  coefficients  of  the  Py 
configurations.  In  order  to  divide  by  (E  -  Hqq)  in  Eq.(13)  and  to  prevent 


the  mixture  of  QV  into  the  Py  that  is  determined  by  diagonalizing  the  H 
matrix  In  Eq.(IO),  It  is  necessary  to  first  solve  for  the  Qf  roots  by 
diagonalizing  Hqq.  The  low  energy  Qy  roots  are  then  projected  out  of  the  Hqq 
portion  of  the  H  matrix  in  Eq.(lO).  Py  is  then  determined  by  diagonalizing 
the  transformed  H  matrix.  The  free  electron  is  represented  by  a  Rydberg 
orbital  with  a  high  principal  quantum  number.  Successive  Py’s  are  determined 
for  a  series  of  increasing  principal  quantum  numbers.’®  The  widths  obtained 
by  this  procedure  are  then  extrapolated  to  the  continuum  to  yield  the 
capture  width. 

These  techniques  have  been  tested  on  the  NO  molecule  and  compared  to 
experimentally  derived**  interaction  matrix  elements  for  the  B  and  L*n 
repulsive  states  of  NO.  These  states  are  known*®  to  play  an  important  role 
in  the  DR  of  the  ground  state  of  NO'*'.  Using  the  approach  outlined  above  with 
a  double  zeta  plus  polarization  Gaussian  basis  set  and  valence  CI*s  of  1335 
terms  for  the  entrance  channel  and  191  terms  for  the  exit  channel  wave 
functions,  comparison  can  be  made  to  the  experimental  results.  The  lowest  4 
roots  in  the  Q  space  have  been  transformed  out  before  solving  for  the 
entrance  channel  wave  functions.  The  experimentally  derived  width  matrix 
elements  have  been  reported  only  for  the  lowest  (n-3,^)  NO  *n  Rydberg  states 
and  not  for  the  high  (n«7,8)  Rydberg  states  used  here  for  extrapolation  to 
the  continuum.  For  the  B*n  repulsive  state  Hamiltonian  matrix  element  with 
the  3pii  C*n  Rydberg  state,  the  calculated  and  experimental®  results  are  1229 
and  1382.6  cm”*,  while  for  the  <4ptt  K*n  Rydberg  state  the  results  are  71^.2 
and  803.9  cm”*,  respectively.  For  the  L*n  repulsive  state,  the  calculated 
and  experimental  matrix  elements  with  the  C  state  are  5^6.7  and  549.0  cm”*, 
respectively.  The  average  difference  of  only  about  8%  between  theory  and 
experiment  for  the  lower  two  levels  is  quite  encouraging,  considering  that 
only  small  valence  space  Cl  wave  functions  have  been  used.  Since  the  matrix 


element  must  be  squared  in  order  to  calculate  the  width,  this  approach  can 
be  expected  to  yield  widths  for  electron  capture  to  an  average  accuracy  of 
17%.  Note  that  this  comparison  of  theory  and  experiment  is  only  for  the 
n-3»^  Rydberg  states  for  which  there  are  experimental  data.  Since  the  n-7,8 
Rydberg  states  are  used  in  the  width  calculations,  an  accuracy  of  better 
than  17%  would  be  expected.  This  approach  has  also  been  used  to  calculate 
widths  for  the  DR  of  0^  along  the  I'llg,  I’Rg,  1*^^,  and  1 ‘Ay  states  of 

Og  and  is  described  in  a  reprint  and  preprint  in  the  Appendix.  The  widths 
range  from  0.8leV  for  ’Ey  to  0.0022eV  for  ’llg. 


IV.  Cross  Sections 


Under  previous  AFOSR  support  the  DR  cross  sections  have  been  calculated 
using  the  quantum  defect  theory  expression**  for  DR  from  ion  vibrational 
level  v: 


agR-(itr/2k|)[Jl^v/(1  +  Iv'5v’)^3  (15) 

where  r  is  the  ratio  of  multiplicities  of  the  neutral  and  ion  states  and  kg 
is  the  free  electron  wave  number.  The  matrix  elements  are  given  by 

Cv"(**/2)[<Xv|rl''2(R)  jXc|>]2  (16) 

where  Xy  and  X^j  are  vibrational  wave  functions  for  the  ion  and  dissociative 
states,  respectively.  r(R)  is  the  width,  discussed  earlier,  that  connects 
the  electronic  continuum  of  the  free  electron  and  ion  to  the  repulsive 


dissociating  resonance  state. 


The  denominator  in  (15)  allows  for  the  reduction  of  the  DR  cross  section 


due  to  autoionization,  i.e.  the  emission  of  the  captured  electron.  The 
summation  over  v'  in  (15)  includes  all  v'  that  are  energetically  accessible 
at  a  particular  electron  energy  and  leads  to  abrupt  decreases  in  the  cross 
section  as  new  channels  open. 

A  computer  program  has  been  written  under  AFOSR  support  which  uses 
expression  (15)  for  the  calculation  of  DR  cross  sections  and  rates.  The 
program  calculates  all  the  required  matrix  elements  using  numerically 
determined  vibrational  wave  functions.  The  calculation  of  the  cross  sections 
and  rates  for  DR  of  Oj  leading  to  0(*S)  has  been  completed  and  a  reprint 
summarizing  the  results  is  Included  in  the  Appendix.  The  calculated  DR  cross 
sections  for  the  lowest  ten  vibrational  levels  of  the  ion  ground  state  are 
shown  in  Fig.  1.  Fig.  la  shows  the  ion  vibrational  wave  functions,  the  ion 
RKR  potential  curve®,  and  the  calculated  * dissociative  route  which  is  the 
only  route  for  generating  0(*S)  from  the  lowest  ion  vibrational  levels. 

Flg.s  1c-k  show  the  large  cross  section  windows.  The  windows  result  from  the 
overlap  of  the  nodes  in  the  ion  vibrational  wave  functions  with  the 
continuum  vibrational  wave  function  in  the  repulsive  state**.  The  number  of 
windows  can  be  predicted  by  simply  counting  the  number  of  nodes  in  the  ion 
vibrational  wave  function  which  are  at  the  small  internuclear  distance  side 
of  the  Intersection  of  the  repulsive  curve  with  the  ion  vibrational  level. 
The  cross  sections  also  show  smaller  abrupt  breaks  that  arise  from 
autolonizatlon.  As  the  electron  energy  is  increased  new  ion  vibrational 
levels  become  accessible  and  can  be  populated  if  autoionization  occurs 
before  dissociation,  leading  to  a  sudden  decrease  in  the  DR  cross  section 
at  the  vibrational  level  thresholds.  The  large  peak  in  the  v«0  cross  section 
(Fig.  1b)  arises  from  the  Franck  Condon  overlap  of  the  continuum  vibrational 
wave  function  of  the  ‘ Ej,  state  with  the  peak  of  the  Gaussian  like  ion  wave 


1 1 


function  for  v-0.  The  calculated  rates  for  DR  are  discussed  further  in  the 
Appendix.  Comparison  of  these  rates  with  the  calculated  and  atmospherically 
derived  rates  found  by  other  researchers  shows  that  the  atmospheric  models 
are  observing  Oj  in  a  non-Boltzmann  vibrational  distribution. 

The  results  for  generation  of  0(‘D)  are  given  in  a  preprint  in  the 
Appendix.  DR  leading  to  0(*D)  from  the  lowest  two  vibrational  levels  is 
dominated  by  the  and  states.  For  the  v=2  level,  the  contribution  of 
is  of  similar  importance  to  that  of  the  other  two  states.  The  total  rate 
from  v-0  is  in  agreement  with  the  rate  derived  from  atmospheric  models.  In 
contrast  to  the  situation  for  0(*S)  the  rate  decreases  with  increasing 
vibrational  excitation  with  the  rate  for  v=2  being  about  one  half  the  rate 
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Fig.  la  The  calculated  potential  curve 
for  the  *Zu  dissociative  route  and  the 
RKR  ion  ground  state  potential  curve 
and  vibrational  wave  functions. 


Fig.  1b  DR  cross  section  for  the  v-0 
vibrational  level  of  the  ion  ground 
state. 


Fig.  1c  DR  cross  section  for  the  v-1 
vibrational  level  of  the  ion  ground 
state. 


Fig.  Id  DR  cross  section  from  v-2. 
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TIm  niMtlnn  that  the  dteMtciative  recombination  (DR)  of  Ot 
nrUh  elactiM  conM  he  ae  important  process  in  the  Earth’s 
■ppcr  atmeapheie  first  appeared  in  the  iiterature  over  55  years 
ato*.  In  1947,  Bates  and  Massey’  pointed  out  that  DR  of  Oj  was 
the  oaly  precses  that  coaid  explain  the  observed  eiectron  recombi- 
nadea  rataa  la  the  E  and  FI  regions  of  the  ionosphere.  In  1954, 
NleaM*  psopoaed  that  DR  of  O2  was  the  source  of  the  ionospheric 
rmlsslnn  at  5,577  A  In  the  Earth’s  airglow.  In  the  intervening  years, 
the  DR  of  leading  to  O  la  the  excited  'S  state  has  been  the 
sa^lect  af  several  reviews  and  away  papers^  and  has  been  termed 
the  ’rinsrical**  taaoapheric  source  of  the  well-known  green  line 
oadaioa.  Neveithclcas,  the  rate  coefficient  for  DR  ieading  to  0('S) 
haa  aevar  heaa  rcliahly  deterasiaed  in  the  laboratory.  Here,  we 
rapert  ths  first  theoretical  calcuiationa  of  the  rate  coefficient  for 
a  wide  range  af  temperatures. 

The  DR  of  O]  with  an  electron  leading  to  excited  atomic 
oxygen  in  the  'S  state  it  described  by 

0j+e--v0+0(‘S)  (1) 

For  the  low  vibrational  levels  of  the  ion  (u<  10),  the  products 
of  reaction  (I)  can  be  described  by  a  single  potential  energy 
curve  of  symmetry’.  In  agreement  with  earlier  theoretical 
studies’,  analyses  of  the  Atmosphere  Explorer  E  satellite 
measurements  have  shown  a  correlation  between  the  quantum 
yield  of  0('S)  from  reaction  (1)  and  the  extent  of  vibrational 
excitation  of  the  ion  although  the  vibrational  distribution  was 
unknown*.  Dissociation  on  a  surface  leading  to  0('D)  +  0('S) 
has  been  supported  bv  interferometer  emission-line  profile 
measurements  at  3,577  A  taken  by  the  Dynamics  Explore  Satel¬ 
lite*.  The  modelling  of  this  process  requires  knowledge  of  the 
dependence  of  the  DR  rate  coefficients  for  production  of  ground 
and  excited  atoms  on  ion  vibrational  excitation  and  electron 
temperature.  Except  for  the  total  DR  rate  coefficient  to  all 
chaMels'*''’,  reliable  rate  coefficients  and  quantum  yields 
leading  to  specific  excited  states  of  O  from  individual  vibrational 
levels  have  not  been  previously  available. 

The  mechanism'*  for  DR  is  illustrated  in  Fig.  1.  The  ion  in 
vibrational  level  v  captures  an  electron  of  energy  e  into  a  neutral 
repulsive  state  which  can  dissociate  before  the  captured  electron 
is  emitted.  The  probability  of  capture  into  vibrational  level  v  is 
proportional  to  a  Franck-Condon  (FC)  factor  between  a  bound 
ion  vibrational  wave  function  and  a  continuum  wave  function 
which  has  a  small  internuclear  distance  (R)  turning  point  near 
the  tip  of  the  vertical  arrow  in  Fig.  1.  As  e  increases,  the  turning 
point  of  the  continuum  wave  function  moves  to  smaller  R,  and 
the  PC  overlap  samples  the  smaller  R  region  of  the  bound  ion 
wave  function. 

The  orbitals  for  the  large-scale  wave  functions  reported  here 
for  the  state  have  been  determined  in  complete  active  space 
self-consistent  field  calculations  (CASSCF)'^  over  a  nuclear- 
centred  contracted  gaussian  basis  set  of  6s,3p,2d,l/  size.  The 
final  configuration  interaction  (Cl)  wave  functions  were  expan¬ 
ded  in  terms  of  the  CASSCF  orbitals  and  generated  from  a 
ten-term  reference  set  which  contained  all  the  configurations 
needed  to  describe  the  dissociation  of  to  a  'D  and  a  'S 
oxygen  atom.  The  Cl  involved  all  single  and  double  excitations 
to  the  (hll  virtual  space  and  yielded  a  139,946-term  wave  func¬ 
tion.  The  SWEDEN  group  of  programs  (written  by  C  W.  Baus- 
chiicher,  P.  E.  M.  Siegbahn,  B.  Roos,  P.  Taylor  and  J  Almhof) 
and  the  direct  Cl  method  of  Siegbahn  were  used  to  determine 
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Fig.  1  Potential  energy  curves  for  the  ground  state  of  O2  (from 

ref.  18)  and  for  the  dissociative  state  calculated  here.  The  ion 

vibrational  wave  functions  for  the  lowest  five  levels  are  each  plotted 
on  the  same  scale. 

the  large-scale  wave  functions'*.  The  final  energies,  corrected 
for  the  missing  quadruple  excitations”,  were  determined  at  58 
points  in  the  range  1.85-8.0  a,,  (Bohr).  The  resulting  potential 
curve  for  is  shown  in  the  region  near  the  ion  in  Fig.  1.  The 
ground-state  ion  potential  curve  is  the  RKR  (Rydberg,  KJein, 
Rees)  curve"  positioned  at  the  experimental  ionization  potential 
above  the  calculated  ground  state  and  shifted  by  0.0180  do  to 
larger  R,  to  compensate  for  the  difference  between  the  calculated 
and  experimental  equilibrium  R. 

The  cross-section,  trfu),  for  DR  through  ion  vibrational  level 
t),  was  calculated  using  the  expression  derived  by  Giusti'*  in 
combination  with  numerically  determined  vibrational  wave 
functions  obtained  by  solving  the  one-dimensional  nuclear 
Schrodinger  equation  in  the  '£„  ^nd  ion  potentials.  The  cross- 
section  is  approximately  directly  proportional  to  an  electronic 
width,  r,  which  couples  the  electronic  continuum  of  the  free 
electron  and  ion  to  the  repulsive  dissociating  resonance  state. 
The  width,  a  matrix  element  of  the  hamiltonian  operator,  is 
calculated  using  Fermi’s  Golden  Rule.  High  principal  quantum 
number  Rydberg  orbitals  are  determined  in  improved  virtual 
orbital’"  calculations  using  an  extended  Rydberg  basis  set  and 
are  used  to  represent  the  free  electron.  The  calculation  of  F  is 
then  a  bound-state  problem.  The  multi-configuration  wave  func¬ 
tion  for  the  ion  plus  high  Rydberg  orbital  (or  free  electron)  is 
optimized  in  the  field  of  an  optical  potential  in  a  procedure 
based  on  the  use  of  Feshbach’’  projection  operators  (S.L.G., 
manuscript  in  preparation).  F  is  determined  for  several  values 
of  the  principal  quantum  number,  n,  of  the  Rydberg  state  and 
a  threshold  value  is  determined  by  extrapolating  to  n-*oo.  The 
resulting  F  calculated  for  the  '£*  state  of  O2  is  0.29  eV.  The  rate 
coefficient,  a,  is  then  obtained  by  averaging  the  cross-sections 
over  a  Maxwellian  distribution  of  electron  energies. 

The  calculated  rate  coefficients  are  shown  in  Fig.  2.  The  values 
of  a  in  the  electron  temperature  ( T,)  range  200  <  T,  <  400  K 
can  be  closely  represented  by  a(t)  =  0)  =  3.2(+1.6,  -l.l)x 
10""’x(r/300)-"'‘,  a(o  =  l)  =  7.3(-l-1.9,  -1.7)xl0-’x 

(r,/300)-"*',  a(u  =  2)  =  2.5(+0.3,  -0.4)  x  10-*x (r./300)-'’‘', 
a(i;  =  3)  =  8,6(  +  5.6,  -4.3)  x  lO’*  x  (r,/300)-'’“  and  a(i!  =  4)  = 
7.2(  +  l  4,  -2.3)xlO-’x(r./300)-"".  The  values  for  a  are  for 
equal  ion  rotational  and  electron  temperatures.  The  precision 
estimates  are  the  average  uncertainties  that  result  for  an  esti¬ 
mated  ±15%  uncertainty  in  the  electronic  width  and  a  ±0.1  eV 
uncertainty  in  the  position  of  the  dissociative  '£*  route  relative 
to  the  ion.  The  u  =  0  level  has  the  lowest  a  and  is  a  factor  of 
78  lower  than  the  highest  or  due  to  u  =  2.  Clearly,  for  the  T, 
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Fig.  1  Hie  calculated  rate  coiuunt  for  diieociative  recombination 
along  the  'IJ  diaaociative  route  for  each  of  the  lowest  five  vibra¬ 
tional  levels  of  the  ion  and  for  the  case  where  the  vibrational  and 
electron  temperatures  are  identical.  In  each  case,  the  ion  rotational 
temperature  is  equal  to  the  electron  temperature. 

range  shown  in  Fig.  2,  the  rate  coefficient  from  p  =  2  is  an  upper 
bound  to  the  coefficient  at  any  vibrational  temperature,  Tvib- 
The  relative  magnitude  of  these  values  can  be  understood  by 
inspection  of  Fig.  t.  At  tow  electron  energies,  the  FC  factor 
between  the  continuum  vibrational  wave  function  of 
ion  vibrational  level  v  increases  with  increasing  o  for  v  <  2.  For 
p  -  0,  the  FC  factor  increases  with  electron  energy  leading  to  a 
rising  above  790  K.  The  increase  is  due  to  the  peak  in  the 
cross-section  at  1.9  eV  which  is  a  'reflection'  of  the  shape  of  the 
p-0  vibrational  wave  function.  Bardsley^^  has  shown  that  a 
variation  of  the  rate  coefficient  as  77*”  arises  if  the  FC  factors 
are  constant  with  electron  energy.  For  v  =  0, 1 , 2, 4  in  the  calcula¬ 
tions  reported  here,  the  FC  factors  increase  as  the  electron 
energy  increases  from  zero  and  the  rate  coefficients  have  an 
exponent  which  is  more  positive  than  -O.S. 

Using  calculated  rates  for  each  of  the  lowest  10  vibrational 
levels  of  the  ion,  values  of  a  for  the  equilibrium  case  in  which 
the  electron  and  vibrational  temperatures  are  identical  have 
been  determined  and  are  shown  in  Fig.  2.  For  T  <  300  K,  the 
rate  coefficient  is  indistinguishable  in  Fig.  2  from  the  coefficient 
for  o«0.  For  r>300K,  the  a  rises  due  to  the  increasing 
contribution  of  the  p>0  rate  coefficients  to  the  total  rate 
coefficient  and  to  the  increase  in  the  v  =  0  rate  coefficient  with 
increasing  T,.  For  800<  T<  1,S00  K,  the  rate  coefficient  is 
closely  represented  by  a  =8.1  x  10"'®x(r/l,150)'^’.  It  is  inter¬ 
esting  to  note  that  atmospheric  models*"''*'’  have  all  used  a 
negative  exponent  for  the  temperature  dependence  of  the  rate 


coefficient  near  1,000  K.  The  positive  exponent  obtained  here 
for  Te  =  near  1,000  K  is  confirmation  that  the  atmospheric 
models  must  be  describing  conditions  for  which  T,  #  Tvii,.  The 
most  recent  atmospheric  model*  reports  the  rate  coefficient  to 
be  (0.83-2.0)  X  10“*  cm*  s"'  for  T,  in  the  range  680  iC-840  K. 
The  results  obtained  here  indicate  that  the  derived  atmospheric 
coefficients  must  pertain  to  O]  which  cannot  be  characterized 
by  a  Boltzmann  temperature  but  which  is  instead  significantly 
populated  in  the  u  >  0  levels. 

The  above  results  indicate  that  a  value  for  the  quantum  yield 
of  0(’S)  is  meaningful  only  if  the  ion  vibrationd  distribution 
is  specified.  Taking  the  ratio  of  the  DR  rate  coefficient  reported 
here  for  p  =  0  to  the  laboratory  value'®  of  the  total  rate  coefficient 
at  T,  =  T,ib  =  300  K  gives  a  quantum  yield  of  0.0016  (+0.0009, 
-O.cibos)  for  0('S)  from  p  =  0  and  should  be  compared  to  the 
original  widely  used  quantum  yield  of  0.1*'.  The  quantum  yield 
reported  here  supports  the  more  recent  view**  that  in  the  earlier 
experiment*'  the  ions  may  not  have  been  vibrationally  relaxed. 
However,  analysis  of  the  experiments**'**  is  difficult  due  to  the 
recent  measurement*'  of  the  rapid  relaxation  of  vibrationally 
excited  O2  in  an  argon  bufler.  The  calculation  of  the  quantum 
yields  from  each  of  the  higher  vibrational  levels  is  currently  in 
progress. 

I  thank  Dr  Richard  1..  Jafle  for  valuable  discussions  concern¬ 
ing  molecular  ion  rotational  considerations.  This  research  is 
supported  at  ISR  by  NASA  Ames  Cooperative  Agreement  and 
at  CFA  by  the  NSF  and  by  AFOSR 

Note  added  in  proof :  In  agreement  with  the  conclusions  reported 
here,  a  new  model  of  the  5,577  A  emission  recorded  by  the 
Dynamics  Explorer  satellite  shows  that  the  p  =  1  and  2  ion 
vibrational  levels  are  the  dominant  contributors  to  0('S)  pro¬ 
duction*®. 
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ABSTRACT-The  results  of  large  scale  ab  initio  calculations  of  the  rates  for 
production  of  0(*D)  by  dissociative  recombination  of  Oj  are  presented  for 
electron  temperatures  in  the  range  IOO-3OOOK.  A  state  is  the  dominant 
dissociative  route  from  v»0  and  a  *1^  state  is  the  most  important  route  from 
v-1  and  v-2.  The  calculated  total  rate  for  0(*D)  production  from  v-0  is 
2.21  (+0.21  ,-0.24)x10“'x(Tg/300)“* ■**  near  room  temperature.  The  v»1  and  v-2 
rates  are  about  ITU  and  ^1%  smaller  respectively,  than  the  v-0  rate  at  300K. 


I.  Introduction 

For  the  molecular  ion,  AB'*',  dissociative  recombination  (DR)  with  an  electron 
is  described  by 

AB""  +  e"  -»  A  +  B  (1  ) 

where  A  or  B  are  ground  or  excited  state  atoms  or  molecular  fragments.  The 
study  of  the  process  of  DR  of  molecular  ions  with  electrons  began  with  the 
suggestion  of  Bates  and  Massey  (19^7)  that  DR  could  be  an  important  sink  for 
electrons  in  the  ionosphere.  It  is  now  recognized  that  DR  can  be  a  rapid 
process  in  a  wide  variety  of  plasma  environments. 

Alex  Dalgarno  has  made  many  significant  contributions  to  the  study  of  DR. 
These  have  ranged  from  fundamental  studies  of  the  DR  mechanism  to  applications 
in  interstellar  and  atmospheric  chemistry.  The  interstellar  studies  (Dalgarno, 
1976,  1982,  1983,  1985,  1986;  Dalgarno  et  al.,  1973;  Black  and  Dalgarno,  1977: 
Dalgarno  and  Black,  1976)  have  explored  the  DR  of  a  wide  variety  of  molecules, 
ranging  in  size  from  to  Several  of  the  atmospheric  studies  of  DR 

(Dalgarno,  1958,  1983)  have  focused  on  the  role  of  DR  in  the  production  of  the 
red  line  of  atomic  0  (Dalgarno,  1958;  Dalgarno  and  Walker,  1964;  Abreu  et  al., 
1986).  Additional  contributions  relevant  to  the  theory  of  DR  have  encompassed 
the  calculation  of  the  energies  (Dalgarno  and  Drake,  1971)  and  widths 
(Bransden  and  Dalgarno,  1952,  1953,  1953a)  of  autoionizing  states  and  the 
definition  of  diabatic  states  (Dalgarno,  1980).  An  analysis  of  the  mechanism 
of  DR  (Bates  and  Dalgarno,  1962),  based  upon  the  earlier  work  of  Bates  (1950, 
1950a),  provided  a  number  of  important  insights  into  DR.  Perhaps  most  relevant 
to  the  work  presented  below  was  the  prediction  that  the  electron  temperature 
dependence  of  the  DR  rate  coefficient  is  Tg®**  unless  the  curve  crossing  of 
the  dissociative  and  molecular  ion  states  leads  to  a  rapid  variation  of  the 
Franck  Condon  factors  with  electron  energy. 

One  of  the  examples  presented  by  Bates  and  Dalgarno  (1962)  was  the  DR  of  0^ 
with  an  electron,  Oj  +  e~  -»  0  +  0.  It  is  now  well  recognized  that  DR  of  Oj  is 
an  important  source  of  0(‘D)  in  the  aurora  (Dalgarno  and  Reid,  1969;  Bates, 
1982;  Rees  and  Roble,  1986)  and  airglow  (Schaeffer  et  al.,  1972;  Bates,  1982; 
Sharp,  1986). 


This  paper  reports  the  first  large  scale  ab  Initio  calculations  of  the  rate 
coefficient  for  generation  of  0(*D)  atoms  from  DR  of  0^  in  one  of  the  lowest 
three  vibrational  levels.  The  next  section  has  a  discussion  of  the  techniques 
used  to  accurately  determine  the  potential  curves  for  the  important 
dissociative  routes  of  Oj.  The  method  used  for  calculating  the  cross  sections 
and  rates  is  discussed  in  Section  III.  The  approach  used  for  determining 
autoionization  widths  is  summarized  in  section  IV  and  the  resulting  rates  for 
0(‘D)  production  are  presented  in  Section  V. 


II.  Potential  Curves 

In  a  seminal  paper,  Bates  (1950a)  pointed  out  that  DR  rates  of  significant 
magnitude  at  low  electron  energies  require  that  a  repulsive  potential  curve, 
describing  the  motion  of  the  products  in  (1),  crosses  the  ion  curve  within  the 
turning  points  of  the  ion  vibrational  level  undergoing  recombination.  Since 
the  time  for  dissociation  can  often  be  rapid  compared  to  the  time  for  emission 
of  the  captured  electron,  the  rate  of  DR  can  indeed  be  fast.  The  process  is 
illustrated  in  Fig.  1  where  an  ion  in  the  v«1  vibrational  level  undergoes  DR 
by  capture  of  an  electron  of  energy  e.  For  the  DR  of  O2,  inspection  of  the 
relevant  potential  curves  (Guberman,  1983)  shows  that  for  the  lowest  three 
vibrational  levels  of  the  ion,  only  five  repulsive  states  can  yield  0(*D).  The 
states  of  interest  are  the  lowest  states  of  ’^u>  *^u>  *^g'  ’'^g 

symmetries.  The  latter  two  states  do  not  dissociate  adiabatically  to  *D  atoms 
but  have  avoided  crossings  with  higher  states  that  lead  to  *D  atoms.  However 
the  two  n  states  have  exceedingly  small  calculated  widths,  .0030eV  for  ‘Ilg  and 
.0022eV  for  ‘llg  (Guberman,  1987a)  at  R»2.28l9ao.  These  widths  are  about  two 
orders  of  magnitude  smaller  than  the  smallest  width  calculated  (see  Section 
IV.)  for  the  other  three  states.  Since  these  small  widths  will  lead  to  very 
small  DR  rates,  these  states  are  neglected  in  the  calculation  of  the  total  DR 
rate  for  production  of  ‘D  atoms.  The  calculated  potential  curves  for  the 
remaining  three  states  are  shown  in  Fig.  1 .  The  potential  curves  have  been 
calculated  using  contracted  Gaussian  basis  sets  of  [63, 3p,2d, 1 f]  size.  The 
basis  sets  have  been  chosen  to  exclude  Rydberg  character  and  the  resulting 
potential  curves  describe  the  diabatic  valence  states  needed  for  DR.  The 
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orbitals  are  determined  separately  for  each  state  In  Complete  Active  Space 
Self-Consistent  Field  (CASSCF)  calculations  (Roos  et  al.,  1980)  and  the 
Configuration  Interaction  (Cl)  wave  functions  are  generated  by  taking  all 
single  and  double  excitations  from  a  multireference  set  of  configurations  to 
the  full  virtual  orbital  space.  The  direct  Cl  method  (Siegbahn,  1980)  and  the 
MOLECULE  (Almlof)-  SWEDEN  (Siegbahn  et  al.)  programs  were  used  for  the 
calculation  of  the  potential  curves.  The  Cl  energies  have  been  corrected  for 
the  missing  quadruple  excitations  (Langhoff  and  Davidson,  197^).  The  final  Cl 
wave  functions  for  the  *1^,  ‘Ay,  and  ‘Zy  states  had  228,036,  185,^00  and 
139,946  terms  respectively.  Further  details  of  the  calculations  will  be 
reported  separately  (Guberman,  1987b). 

The  *Zy  state  shown  in  Fig.  1  Intersects  the  ion  at  the  small  R  turning 
point  of  the  v«0  level  and  is  the  dlabatic  continuation  of  the  repulsive  wall 
of  the  upper  state  of  the  well  known  Schumann-Runge  bands.  Dissociation  along 
*Zu  after  capture  by  a  zero  energy  electron  from  v-0  leads  to  a  ground  state 
atom  and  a  *D  atom,  each  with  2.5  eV  kinetic  energy.  Dissociation  along  the 
‘Ay  route  from  v-0  leads  to  two  ‘D  atoms,  each  with  1.5eV  kinetic  energy. 
Dissociation  along  ‘Zy  leads  to  a  *S  and  a  *D  atom,  each  with  0.4  eV  kinetic 
energy  from  v-0. 


111.  Methods  for  Cross  Sections  and  Rates 

The  cross  sections,  Oy,  for  direct  DR  (Bardsley,  1968)  from  ion  vibrational 
level  V  have  been  calculated  using  the  expression  derived  by  Giusti  (1980) 
from  quantum  defect  theory  for  the  case  of  a  single  dissociative  state  and 
many  open  ion  vibrational  levels,  v’, 

Ov-(TTr/2k|)[4^v/(1+Iy,Cv')2]  (2) 


where  kg  -  2me ,  kg  is  the  wave  number  of  the  incident  electron, 
electron  mass,  and  r  the  ratio  of  multiplicities  of  the  neutral 
states.  The  matrix  elements  Cy  given  by 


m  is  the 
and  ion 
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Cv-(’f/2)[<Xv|r1/2(R)  |Xd>]2 


(3) 


where  Xy  and  X^j  are  vibrational  wave  functions  for  the  ion  and  dissociative 
states,  respectively.  The  integration  in  (3)  is  over  both  R  and  the  electronic 
coordinates.  The  vibrational  wave  functions  are  determined  by  the  Nuraerov 
method  on  a  O.OOIao  grid  between  0.0  and  8.0ao.  The  integration  over  R  in  (3) 
is  done  by  Simpson's  rule.  The  width,  r{R),  connects  the  electronic  continuum 
of  the  free  electron  and  ion  to  the  repulsive  dissociating  resonance  state  and 
is  given  by  the  Golden  Rule, 

r(R)-2Tr<{*ion(x,R)^£.(x,R)}|H|*d(x,R)>2.  (U) 

In  (4)  H  is  the  Hamiltonian,  the  integration  is  over  the  electronic 
coordinates  represented  by  x,  and  the  wave  functions  from  left  to  right  in  (4) 
are  for  the  ion  core,  the  free  electron,  and  the  dissociative  neutral  state. 
The  antisymmetrized  product  of  the  ion  core  and  free  electron  wave  functions 
is  denoted  by  {}. 

The  indirect  DR  process  (Bardsley,  1968a;  Giusti,  1980;  Giusti-Suzor  et  al., 
1983)  in  which  recombination  occurs  into  a  vibrationally  excited  Rydberg  state 
which  is  predissociated  by  the  repulsive  routes  discussed  here,  has  not  been 
Included  in  the  cross  section  calculations.  The  indirect  process  is  not 
expected  to  be  a  significant  contributor  to  0(*D)  production. 

Rate  constants  have  been  calculated  by  averaging  the  calculated  cross 
sections  over  a  Maxwellian  distribution  of  electron  energies. 

IV.  Electronic  Widths 

The  widths  given  in  Eq.  (4)  have  been  calculated  by  using  high  Rydberg 
orbitals  to  represent  the  free  electron  (Guberman,  1986,  1987a).  Since  all  the 
orbital  amplitudes  in  (4),  except  for  that  of  the  free  electron,  fall  off 
rapidly  with  increasing  x,  only  the  free  electron  amplitude  near  the  nuclei  is 
needed  for  the  width  calculation.  Except  for  a  normalization  factor,  high 
Rydberg  orbitals  have  a  shape  similar  to  free  electron  orbitals  near  the 
nuclei.  A  density  of  states  is  then  Inserted  into  (4)  for  the  normalization 
(Guberman,  1987a).  The  Rydberg  orbitals  are  calculated  using  the  Improved 
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Virtual  Orbital  (IVO)  method  (Hunt  et  al.,  1969)  in  a  basis  set  of  18  diffuse 
Gaussians  positioned  midway  between  the  nuclei  (Guberman,  I987a).  Each  Rydberg 
orbital  is  an  eigenfunction  of  an  IVO  Hamiltonian  which  is  unique  for  each 
molecular  symmetry. 

The  wave  function  of  the  ion  plus  a  free  electron  in  (4)  is  represented  by 
antisymmetrized  products  of  a  valence  Cl  description  of  the  ground  state  of  0^ 
with  the  Rydberg  orbital  and  with  the  valence  virtual  orbitals  of  the  same 
symmetry  as  the  Rydberg  orbital.  For  the  width  calculations,  the  dissociative 
states  were  represented  by  full  valence  Cl  wave  functions  in  which  excitations 
out  of  the  2o  orbitals  were  included.  For  this  procedure  leads  to  a  ^8 
term  Cl  wave  function  for  the  dissociative  state  and  a  954  term  Cl  wave 
function  for  the  entrance  channel.  For  ‘Ay  and  entrance  and  exit 

channel  wave  functions  have  42  and  630  terms  respectively.  The  size  of  the 
entrance  channel  wave  function  is  Increased  further  by  allowing  the  entrance 
channel  to  mix  in  higher  nonphysical  routes  from  the  exit  channel.  A  Feshbach 
(1958)  projection  operator  procedure  has  been  used  to  implement  this  procedure 
(Guberman,  1986,  1987a). 

The  entire  procedure  is  repeated  for  several  values  of  the  Rydberg  orbital 
principal  quantum  number  with  the  highest  accurately  determined  orbital  having 
n-8.  The  widths  calculated  from  each  value  of  n  are  then  extrapolated  to  the 
continuum.  A  comparison  to  experimentally  derived  (Galluser  et  al.,  1982) 
predissociation  matrix  elements  for  NO  Indicates  that  the  widths  calculated  by 
this  procedure  deviate  from  experiment  by  about  15)1. 

For  *Eu,  the  calculated  width  is  O.81  eV  near  R-2.08ao  and  decreases  to 
0.54eV  at  1.78ao.  For  ‘A^  and  *Ey,  the  width  is  relatively  constant  with  R, 
having  a  value  of  0.44eV  and  0.29eV  respectively  at  R-2.28ao .varying  by  only 
about  2%  over  the  region  of  significant  nuclear  overlap  with  the  v-0,1,2  ion 
wave  functions.  Since  the  estimated  uncertainty  in  the  calculated  widths  is 
15){,  the  variation  of  the  width  with  R  has  been  ignored  for  'Ay  and  'Ey  but 
has  been  accounted  for  in  expression  (3)  for  ’Ey  (S.  L.  Guberman,  1987b). 

V.  Calculated  DR  Rates 

The  calculated  rate  constants  for  the  three  dissociative  routes  discussed 
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here  are  shown  in  Fig.s  2,  3t  and  for  each  of  the  lowest  three  vibrational 
levels  of  the  nonrotating  ion.  For  200iTeS**00K,  the  DR  rates  are  given  in 
Table  1.  The  uncertainties  in  the  rates  are  average  uncertainties  computed 
from  an  estimated  ±15J  uncertainty  in  the  calculated  widths  and  an  estimated 
±0.1 eV  uncertainty  in  the  position  of  the  potential  curves  relative  to  the  ion 
curve.  The  highest  DR  rates  are  from  *1^.  However,  for  v»0,  the  rate  along  My 
is  more  than  half  the  DR  rate  along  ‘Zy.  Since  leads  to  two  ‘D  atoms,  the 
highest  rates  for  producing  *D  atoms  from  v»0  are  from  ‘Ay.  However,  from  v-1 
and  v-2  (Fig.s  2-4),  the  highest  rates  for  0(*D)  production  are  from  *1^. 

As  pointed  out  by  Bates  and  Dalgarno(1962) ,  the  temperature  dependence  of 
the  direct  DR  process  will  vary  as  Te“**  if  the  amplitude  of  the  ion  wave 
function  does  not  vary  substantially  over  the  range  of  electron  energies  of 
Interest.  Such  a  situation  can  arise  for  any  ion  wave  function  if  the 
dissociative  curve  is  very  steep.  In  general,  an  intersection  of  the 
dissociative  route  with  the  ion  vibrational  level  near  the  peak  of  the  ion 
wave  function  will  lead  to  slowly  varying  nuclear  overlap  terms  and  electron 
temperature  dependencies  near  -0.5,  while  intersection  at  a  point  where  the 
absolute  value  of  the  ion  wave  function  amplitude  is  increasing  (decreasing) 
will  usually  lead  to  temperature  dependencies  more  positive  (negative)  than  - 
0.5.  Note  that  the  denominator  in  (2)  has  a  tempering  effect  upon  these 
conclusions.  Two  such  cases  are  illustrated  in  Cases  a  and  b  of  Fig.  2  of 
Bates  and  Dalgarno  (1962).  The  latter  Case  shows  a  dissociative  route 
intersecting  the  v-0  level  in  the  nonclassical  region  on  the  large  R  side  of 
the  ion  equilibrium  internuclear  distance  and  resembles  the  situation  found 
here  for  ‘Zy. 

For  ‘Zy,  the  amplitudes  of  the  ion  wave  functions  at  low  electron  energies 
Increase  most  rapidly  for  v-0.  The  Increase  for  v-1  is  less  rapid  than  for 
v-0.  For  v-2,  the  intersection  is  near  the  outer  peak  in  the  v-2  wave 
function.  The  calculated  temperature  dependencies  given  in  Table  1  for  *Zy 
show  the  predicted  pattern.  The  most  positive  exponent  is  for  v-0  followed  by 
v-1  and  v-2.  For  the  other  two  dissociative  routes  the  exponents  are  near  0.5 
for  v-0  since  the  intersections  with  the  ion  are  within  the  classical  turning 
points  and  the  amplitude  is  changing  slowly.  For  v-1  and  2  and  *Zy  and  'Ay, 
one  must  take  into  account  the  overlap  of  the  nodes  in  the  ion  wave  functions 
with  the  nodes  in  the  continuum  wave  functions.  The  exponents  are  difficult  to 


predict  in  these  cases  by  visual  inspection  of  Fig.l, 

Also  shown  in  Fig.s  2-^  are  the  rates  for  the  equilibrium  case  where  the 
vibrational  temperature  is  set  equal  to  the  electron  temperature.  For  ’Zu*  the 
equilibrium  case  rates  are  nearly  indistinguishable  from  the  v-0  rates  over 
the  span  of  electron  temperatures  plotted.  For  the  equilibrium  rates  fall 

below  the  v-0  rates  at  high  electron  temperatures.  For  *1^1  the  equilibrium 
rates  increase  dramatically  above  450K,  since  the  rates  for  v-1  and  v-2  are 
substantially  larger  than  for  v-0  (Guberman,  1987). 

It  is  interesting  to  examine  the  ‘D  production  in  each  channel  as  a  function 
of  ion  vibrational  level.  From  v-0,  56%  of  the  *D  atoms  arise  from  the  ‘Ay 
channel  and  ^4$  arise  from  the  *Zu  channel.  For  v-1,  58%  of  the  *D  atoms  are 
in  the  *Zy  channel  and  4$  appear  in  the  *Z^  channel  with  0.52eV  kinetic 
energy.  For  v-2,  45*  of  the  atoms  are  in  the  2.7eV  ’Zy  channel  and  21*  are  in 
the  0.64eV  *1^  channel.  Each  vibrational  level  will  lead  to  a  unique  three 
component  Doppler  line  shape.  From  these  data,  the  vibrational  distribution  of 
the  recombining  ion  could  be  determined  in  an  experiment  that  accurately 
determines  the  Doppler  widths  of  the  radiating  *D  atoms. 

The  total  rate  for  producing  ‘D  atoms  for  Tg  in  the  above  range  from  v-0 
(Table  1)  is  2.21  (+0.21  ,-0.24)  x  10“’  x  (Te/300)"®’‘'‘cmVsec  and  is  in  good 
agreement  with  the  value  originally  reported  by  Zipf  (1970)  of  1 .9  x  10“^ 
cm*/sec.  Note  that  the  total  calculated  rate  decreases  with  increasing  v.  The 
rate  from  v-2  is  about  half  the  rate  from  v-0.  The  variation  of  the  rate  with 
V  contrasts  sharply  to  that  for  producing  0(*S)  (Guberman,  1987)  shown  in  the 
fourth  column  of  Table  1  and  Fig.  4.  The  rotationless  0(*S)  rate  for  v-2  is  a 
factor  of  80  larger  than  the  rate  from  v-0.  The  total  calculated  rate  up  to 
Tg-SOOOK  for  *D  production  is  shown  in  Fig.  5.  Rates  for  both  J-0  and 
rotationally  averaged  rates  in  which  the  rotational  temperature  is  the  same  as 
Tg  are  shown.  Clearly,  accounting  for  the  rotational  temperature  has  only  a 
slight  effect  upon  the  calculated  rates.  The  rotationally  averaged  and  the 
rotationless  results  from  v-0  and  for  the  case  of  equal  vibrational  and 
electron  temperatures  are  each  superimposed  on  the  scale  of  Fig.  5. 

There  have  been  several  aeronoraic  determinations  of  the  *D  production  rate 
by  DR.  The  literature  prior  to  1982  has  been  reviewed  by  M.  R.  Torr  and  D.  G. 
Torr  (1982).  Among  the  deduced  production  rates  are  (2 . 1 ±0. 48)x1 0“’x 
(Tg/300)“-**  (P.  B.  Hays  et  al,  1978;  R.  Link  et  al. ,  1981)  and 
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(1 .96±0.67)x10“’x(Tg/300)“* **  (Abreu  et  al.,  1983).  These  rates  were  obtained 
by  combining  the  total  rate  found  by  Walls  and  Dunn  {197*1)  and  Torr  et  al. 
(1976)  with  a  model  determination  of  the  *D  quantum  yield.  A  recent 
determination  (Abreu  et  al. ,  1986)  which  improves  the  models  by  including 
quenching  of  0(*D)  by  0(*P)  leads  to  a  rate  of  1 .9x10~’x(Tg/300)“' **.  The 
deduced  rates  are  in  excellent  agreement  with  the  v-0  rate  calculated  here  at 
room  temperature.  The  atmospheric  quantum  yields  are  derived  at  ionospheric 
temperatures.  At  Tg-760K,  the  *D  production  rate  calculated  here  is 
1 .41  (+0.1*1, -0.l6)xl0~’  compared  to  the  model  rate  (Abreu  et  al.,  1986)  of  1.14 
X  10“’.  Considering  the  likely  precision  on  the  latter  rate,  the  rates 
reported  here  are  in  agreement  with  the  model  rates.  Note  that  the  precision 
on  the  rates  calculated  here  do  not  exclude  the  possibility  that  the  models 
apply  to  vibrationally  excited  Oj. 

None  of  the  above  models  have  included  ‘D  production  due  to  cascade  from  'S. 
The  DR  contribution  to  cascade  will  be  unimportant  if  the  ion  is  mostly  in  v-0 
but  should  be  considered  if  a  substantial  fraction  of  the  ions  are 
vibrationally  excited. 

VI.  Conclusions 

Three  neutral  dissociative  states  of  *1^,  ‘Ay,  and  * symmetries  have  been 
identified  as  the  major  sources  of  0(‘D)  from  DR  of  the  lowest  three 
vibrational  levels  of  0^.  DR  rates  along  each  of  these  routes  have  been 
calculated  for  electron  temperatures  in  the  range  IOO-3OOOK.  The  total  rate 
for  generating  DR  from  v-0,  2.21  (+0.21  ,-0.24)  x  10“’  x  (Tg/300)“® • '* ‘cm*/sec, 
agrees  with  prior  rates  determined  from  aeronomic  models  and  laboratory 
measurements.  However,  the  variation  of  the  rate  with  v  and  the  precision  on 
the  calculated  rates  leaves  room  for  the  possibility  that  the  atmospheric 
models  and  the  laboratory  experiments  could  be  observing  vibrationally  excited 
ions.  The  total  rate  for  generating  0(*D)  decreases  with  Increasing  v  in 
contrast  to  the  rates  for  generating  0(‘S)  which  increase  markedly  with 
increasing  v. 
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TABLE  1.  Dissociative  recombination  rates  along  three  channels  for  each  of  the  lowest  three  vibrational  levels  of  the 
ion  at  200KST-SJ400K  and  J-0.» 
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Figure  Captions 


1.  Potential  energy  curves  for  DR  leading  to  0(‘D).  The  three  important 
calculated  repulsive  routes  for  DR  of  the  lowest  three  vibrational  levels  are 
shown.  The  bound  potential  for  the  ion  is  taken  from  Krupenie  (1972)  and 
calculated  points  (Guberman,  to  be  published).  The  vibrational  wave  functions 
are  each  plotted  on  the  same  scale  and  are  obtained  from  the  numerical 
solution  of  the  nuclear  Schroedinger  equation. 

2.  Calculated  DR  rates  along  the  *Zu  dissociative  route  leading  to  0(*P)  + 
0(*D)  from  the  lowest  three  vibrational  levels  of  the  ion.  Also  shown  is  the 
rate  for  the  case  in  which  the  vibrational  temperature  is  the  same  as  the 
electron  temperature.  The  rates  are  for  non-rotating  molecules. 

3.  Calculated  DR  rates  along  the  ‘Ay  dissociative  route  leading  to  0(*D)  + 
0(*D)  from  the  lowest  three  vibrational  levels  of  the  ion.  Also  shown  is  the 
rate  for  the  case  in  which  the  vibrational  temperature  is  the  same  as  the 
electron  temperature.  The  rates  are  for  non-rotating  molecules. 

Calculated  DR  rates  along  the  dissociative  route  leading  to  0(*S)  + 
0(*D)  from  the  lowest  three  vibrational  levels  of  the  ion.  Also  shown  is  the 
rate  for  the  case  in  which  the  vibrational  temperature  is  the  same  as  the 
electron  temperature.  The  rates  are  for  non-rotating  molecules. 

5.  Calculated  total  DR  rates  for  production  of  0(*D)  from  the  lowest  three 
vibrational  levels  of  the  ion.  Both  the  J-0  rates  and  rates  for  which  the 
rotational  temperature  is  the  same  as  the  electron  temperature  are  shown. 
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